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Abstract 

We study truncation effects in the SU(3) gauge actions obtained by the Monte Carlo renor- 
malization group method. By measuring the heavy quark potential we find that the truncation 
effects in the actions coarsen the lattice by 40-50% from the original blocked lattice. On the 
other hand, we find that rotational symmetry of the heavy quark potentials is well recovered 
on such coarse lattices, which may indicate that rotational symmetry breaking terms are easily 
cancelled out by adding a short distance operator. We also discuss the possibility of reducing 
truncation effects. 



1 Introduction 



Recently a lot of attention has been devoted to improvements of lattice discretized actions. There 
exist two approaches to improving actions. One is the perturbative improvement program suggested 
by Symanzik|l|] and the other the renormalization group improvement program by Wilson M. Early 
attempts at the perturbative improvement program did not appear practical for Monte Carlo 
simulations j|] . Recently, however, it has been revitalized with the help of tadpole improvement^] 
and actively been investigated by Monte Carlo simulations. 

The renormalization group improvement program is very attractive since it can, in principle, 
give us lattice-artifact-free actions. In practice, however, it is not an easy task to obtain such 
actions since we do not know a practical way to determine them. Recent successful attempts |5| 
approximate the perfect action, which is defined on the renormalized trajectory, with the fixed 
point action ( classically perfect action ) obtained in the limit of (3 — > oo. An advantage of this 
method is that the fixed point action is rather easily obtainable. Although the fixed point action is 
an approximation to the perfect action, this approximation turns out to be rather good in Monte 
Carlo tests. 

Direct attempts^], ||, 11, 12] to obtain actions on the renormalized trajectory have been also 
made, using the Monte Carlo renormalization group (MCRG) methodQ. The MCRG method is 
useful not only for determining improved actions but also for other purposes, like monitoring the 
flow of couplings under a scale transformation. In lattice QCD for instance, the MCRG method 
has been successfully used for determining the coupling shift A/3 of the f3 function to reveal the 
scaling behavior of the theory ||. 

The MCRG approach uses the fact that the blocked trajectory generated by the blocking trans- 
formation reaches the renormalized trajectory after sufficiently many blocking steps and then runs 
along with the renormalized trajectory. Therefore the configurations generated by successive block- 
ing transformations will correspond to an action located nearer and nearer the renormalized tra- 
jectory. Such configurations should have less artifacts than non-blocked configurations since from 
Wilson's renormalization group argument we expect that continuum physics is realized on the 
renormalized trajectory. Actually it has been shown that rotational symmetry of the heavy quark 
potential is well recovered on blocked configurations P]. If we can determine the action represent- 
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ing the blocked configurations, we can directly generate the configurations without blocking. The 



determination of the action was tackled by the canonical demon method [|Tl|, 12 1 which produces as 
output a set of coupling constants. In the canonical demon method, like other determination meth- 
ods, the action to be determined must be truncated to a certain local form, which may cause some 
error. We call the effects of this error truncation effects. Unless the truncated form of the action 
is sufficiently close to the real one, truncation effects may affect long-distance, physical properties 
of the action. In Ref|| several SU(3) gauge actions corresponding to blocked configurations (gen- 
erated by a blocking transformation) were obtained in multi-dimensional coupling space. However 
it is not known whether these actions still preserve the improvements of the blocked configurations 
or not, since these improvements may have been ruined by truncation effects. This point can only 
be examined by Monte Carlo simulations. In this letter, we perform Monte Carlo simulations with 
the actions obtained in Ref||, which we call MCRG improved actions, examine the rotational sym- 
metry of the heavy quark potential, and estimate the truncation effects which appear in a physical 
quantity ( the string tension ). 

2 Monte Carlo Renormalization Group Improved Actions 

We briefly sketch the method used to determine the MCRG actions in Ref||. First, we block 
configurations generated with the standard Wilson action on 32 3 x 64 lattices. The blocking scheme 
employed is Swendsen's scale factor 2 blocking scheme. This blocking scheme was optimized by 
multiplication by a Gaussian random SU(3) matrix so that the blocked trajectory converges to the 
renormalized trajectory quickly [ffO, [9|. Actually the optimal width of the Gaussian distribution at 
[3 ~ 6.0 turned out to be 0, meaning that the Gaussian random SU(3) matrix becomes unity. This 
optimal width was used in the blocking. The blocking was performed twice, starting from 32 3 x 64 
lattices at two (3 values, 6.0 and 6.3, thus resulting in two sets of 8 3 x 16 blocked configurations. 



Next, in order to determine coupling constants we apply the canonical demon method [11, 12 
on the blocked configurations. The canonical demon method introduces several degrees of freedom, 
so-called demons, which are associated with each of the coupling constants to be determined. The 
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action S can be written as 

S = Y^PiSi(U) (1) 

i 

where /3j is a coupling constant and Si(U) is some operator consisting of Wilson loops. Hereafter 
for simplicity, let us assume that the action is characterized by only one operator or by one coupling 
constant: 

S = /3S(U). (2) 

The demon is updated by a microcanonical simulation in the joint system, i.e. the demon and the 
links of a blocked configuration. In the microcanonical simulation, the total energy of the joint 
system, i.e. S(U) plus the demon energy E d , is kept constant. 

In the canonical demon method, in order to avoid a possible finite volume error, a set of well 
uncorrelated blocked configurations is prepared. At a certain stage of the microcanonical simulation, 
we move to a new blocked configuration chosen from the set and the configuration used before is 
discarded. The demon is also moved to the new configuration keeping the value of the demon 
energy at the last update. 

The probability distribution of the demon energy P(E d ) in the simulation is expected to be the 
Boltzmann distribution: 

P(E d ) ~ exp(-/3£ d ) (3) 
where j3 is the coupling constant to be determined. Using eq.(Q), we write 

<E d >= E d eM-PE d )dE d /Z (4) 

where the demon energy E d is restricted in the region, E m % n < E d < E max and Z is the partition 
function. If we take E max = —E m i n = E c , where E c stands for some constant value which we fix 
in the simulation, eq.(H) will be 

< E d >= 1/(3 - E c /tanh{l3E c ). (5) 

Finally, substituting the average value < E d > obtained from the simulation, we determine the 
value of the coupling constant (3 by solving eq. (|5|) numerically. The extension to the multi-coupling 
form of eq.(||) is straightforward. 
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In our study we take MCRG actions obtained in two-dimensional coupling space. The actions 
are written as 

5 = 22e(/3n5^rr(l x 1 Wilson loop)/3 + [3 l2 ^ Tr(l x 2 Wilson loop)/3). (6) 

The values of the couplings (3\i and j3\2 are listed in Table 1. The actions A and B come from RefQ , 
which are obtained from the configurations blocked twice at (3 = 6.00 and (3 = 6.30 respectively. 
We also use an interpolated action (M) located half-way between actions A and B. Fig.l shows the 
locations of these actions in the /3n-/?i2 plane. The ratio of the two couplings, |/?i2//?n|, is bigger 
than that of the Symanzik tree-level improved actions [fi3|1 and Iwasaki action^]. 

3 Heavy Quark Potentials 

We employ 8 3 x 16 lattices which is the same lattice size with the original blocked lattices. We 
calculate static potentials between a heavy quark and antiquark pair from the exponential fall-off 
of Wilson loops. Our calculation is based on 500-700 configurations separated by 100-200 pseudo 
heat-bath sweeps. The smearing technique is used to reduce errors in the extracted potentials. 

By construction, the canonical demon method preserves the average value of the Wilson loops 
included in the action. This means that in our case the average values of 1 x 1 and 1x2 Wilson 
loopsQ, which are associated with j3n and /3i2, should be the same values as the blocked ones within 
statistical uncertainty. Fig. 2 and 3 show comparisons of 1 x T Wilson loops between the original 
blocked lattices and the lattices obtained with the actions A and B. As seen in both figures, good 
agreement is observed for each of 1 x 1 and 1x2 Wilson loops. This check ensures that the canonical 
demon method worked correctly. For other loops, however, the difference increases as T increases. 

In order to compare the potentials we plot the potentials of the action A(B) and original 
blocked lattice on the same figure(Fig.4-5). It is clearly seen that the potentials of the action A(B) 
are very different from the ones of the original blocked lattice, and the slopes of the potentials at 
large distance of the action A(B) are larger, which implies that the lattice spacings are also larger 
than the original ones. We quantify this effect by comparing string tensions. 

Some average values of Wilson loops of the MCRG actions are listed in Table 2. 
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The string tensions are extracted by fitting the potentials to the form 

Q 

V L (r) =m + a L r , (7) 

r 

where m, ctl and c are fitting parameters, and r = R/a is the distance measured in lattice units. 
The fits are performed on all the data including on and off axis potentials. The fit results are 
summarized in Table 3. 

In order to evaluate the truncation effect, let us compare the string tensions with the ones of 
the original blocked lattices. The results are the following: 

1.29(12) MCRG Action A 
°L={ (8) 
0.58(4) Original Blocked Lattice(at p = 6.00) 



0.587(27) MCRG Action B 

(9) 

0.287(26) Original Blocked Lattice(at (3 = 6.30). 



Here we comment on the naive expected string tension on a blocked lattice at /3 = 6.00 and 
6.30. From the literature [DJ, we find aa 2 ~ 0.051 at j3 = 6.00 and aa 2 ~ 0.0200 at = 6.30. Since 
we block twice, the naive expected string tension on the blocked lattice will be aa 2 ~ 0.82(0.32) 
at f3 = 6.00(6.30) respectively. These values are compatible with that of the blocked lattice at 
(3 = 6.30 but not compatible at (3 = 6.00 ( See eqs.(|8|) and (||) ). Probably this mismatch is due to 
the very small correlation length at (3 = 6.00 since the correlation length £ is estimated to be £ ~ 4 
using l/£ 2 ~ aa 2 and this correlation length is too small to preserve the same long range physics 
under the blocking transformation. 

We examine the truncation effects which are seen in the string tension by taking the ratio of 
the string tensions of the MCRG actions and the original blocked lattices (eqs. H[||). 

(o-l)mcrg j 2 - 22 MCRG Action A 
{(Tl) Blocked ) 2.04 MCRG Action B. 

2 We interpolate the value of the string tension for [3 = 6.30. since we do not find any Monte Carlo result at 
P = 6.30. 
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Since ol = aa 2 , the ratio in lattice spacing will be: 

(a) MCRG j 1.5 MCRG action A 
(a) Blocked j 1.4 MCRG action B. 



(10) 



If the truncation effects were negligible, the ratio should be one. However this is far from being 
the case, which indicates that sizeable truncation effects are involved in the MCRG actions. It is 
important to notice here is that the truncation effect increases the lattice spacing, i.e. the lattice 
spacing of the MCRG actions is 40 — 50% bigger than that of the original blocked lattice. This 
increase can be understood in the following way. The MCRG actions contain only the coupling 
constants associated with small Wilson loops (lxl and 1x2 Wilson loops). On the other hand, 
the real action of the blocked configurations can have many coupling constants associated with 
large Wilson loops. Let us assume that all the long range coupling constants are positive. This 
assumption seems to be valid at least for coupling constants up to eight-links Wilson loop operators 
as found in ||, where all those coupling constants are positive. The positive coupling constants 
constrain Wilson loops to be more ordered. Since the large Wilson loops of the MCRG action 
are less ordered average values of the Wilson loops decay faster with distance than that of the 
constrained Wilson loops on the blocked configurations, as seen in Fig. 2-3. Thus the potential of 
the MCRG action, V{R)mcrg-, win be- larger than that of the blocked configurations, V(R)uock : 

V{R)mcrg > V(R) bl ock- (11) 

This situation may become more pronounced at large R as seen in Fig. 4-5. Thus, we obtain a larger 
string tension as extracted from the slope of the potential at large distance. 

In order to examine the rotational symmetry of the potentials we plot all the data of the MCRG 
action A,B and M on the same figure ( Fig. 6 ) by rescaling them using the obtained string tension 
results. The potentials are shifted so that each fitted curve for V{R)/\fa gives the value 2 at 
Ry/a = 2. We see good rotational symmetry except at large R\fa due to the large error bars. The 
recovery of rotational symmetry will be more remarkable if we compare the potentials with that of 



the standard Wilson action [16] where rotational symmetry is severely violated on a coarse lattice. 
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4 Discussion 



The MCRG improvement approach has severe difficulties in determining the effective action. When 
we determine the action, we have to truncate it to a certain local form. On the other hand the 
truncation effects may ruin some benefits from the MCRG improvement. In order to reduce the 
truncation effects, one can tune the blocking transformation in such a way that contributions of long 
range coupling constants disappear quickly. This tuning was successfully done for the fixed point 
action^]. At finite j3, however, a practical way to tune the blocking transformation is not known. 
The brute force approach, i.e. to directly search for an optimal blocking scheme by varying a couple 
of blocking parameters in Monte Carlo simulations, would certainly require a huge computational 
effort which we try to avoid. 

Our MCRG improved actions have been obtained with a blocking scheme optimized by a Gaus- 
sian random SU(3) matrix. Strictly speaking, however, this optimization does not guarantee that 
the truncation effect will be minimized. Originally this optimization was introduced for the study 
of the A/3, and the optimization was done so that the blocked trajectory converges to the renor- 
malized trajectory quickly. This quick convergence is only desirable for the matching method of 
the A/3 analysis. In this sense it is not sure that our MCRG actions are characterized by a local 
form. Actually our study found a big difference in the lattice spacing by the analysis of the string 
tension. Even so good news were found: one of the evidences for improvement, recovery of the 
rotational symmetry, is observed even in the presence of the truncation effects. 

The one-loop level tadpole improvement scheme[Q] is known to work well, in which the action is 
designed to remove the 0(a 2 ) and 0{a s a 2 ) errors by adding two additional operators to the stan- 
dard Wilson action. As far as rotational symmetry is concerned, the tree-level tadpole improvement 
scheme ( Symanzik action + tadpole ) also works well, in which one additional operator (1x2 
Wilson loop) is added to the standard Wilson action. The only difference between the tree-level 
improvement ( Symanzik action ) with and without the tadpole scheme is that the tadpole scheme 
increases the negative coupling contribution of the 1x2 Wilson loop operator, i.e., if we use the 
average plaquette value < plaq >~ 0.4, /?i2//3n is approximately equal to —0.08, instead of —0.05 
for the tree-level improvement scheme without tadpole improvement. For our MCRG actions, we 
find Pu/(3n — —0.1. The Iwasaki action, with which good rotational symmetry can be seen [ 17 1, 
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also has similar behavior, i.e. Piu/Pu — —0.09. Therefore we suggest that rotational symmetry 
can be easily recovered with a short-distance correction (1x2 Wilson loop ) to the standard 
Wilson action and the value of /3i2//?ii should be more negative than the naive tree-level value, 
i.e., 012/Pxi = —0.08 to —0.1 in a region where the lattice spacing a is as large as a ~ 0.3 to 0.5 
fm. 

The MCRG improved actions were determined with the canonical demon method. One could 
try another scheme such as the Schwinger-Dyson (SD) equation method) IS]. Truncation effects 
with the SD equation method might be different, since in that case coupling constants are obtained 
by solving a set of equations whose coefficients are calculated from correlations among arbitrary 
Wilson loops, while the canonical demon method obtains coupling constants so that average values 
of Wilson loops included in the action only are conserved. And even though the canonical demon 
method conserves the average values of such Wilson loops, truncation effects already appear in the 
correlation among them|12]. Thus the SD equation might be more advantageous than the canonical 
demon method since the SD equation uses additional information on the action from the correlation 
among Wilson loops and determines the coupling constants by using many available equations. For 
the SU(3) gauge theory, determination of the renormalized trajectory was attempted with the SD 



equation method[19| and a rough picture of the RT was obtained. It would be interesting to see how 
truncation effects appear in a different way for both the canonical demon and the Schwinger-Dyson 
equation methods. 
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action 


0n 


012 


A 


6.1564 


-0.6241 


M 


7.0712 


-0.7705 


B 


7.986 


-0.9169 



Table 1: Values of coupling constants, (3u and (3±2, used in the study. We label these actions by 
A, M and B as in the table. 



Action 


A 


M 


B 


I x J 


1 2 


1 2 


1 2 


1 


0.409306(75) 


0.475636(77) 


0.538690(60) 


2 


0.139237(62) 0.010488(55) 


0.19660(10) 0.026427(54) 


0.263218(86) 0.056835(92) 


3 


0.045152(65) 0.004230(39) 


0.071291(95) 0.003129(62) 


0.127497(95) 0.01255(10) 


4 


0.014484(53) 


0.036781(82) 0.000404(51) 


0.061726(93) 0.002846(89) 


5 


0.004591(73) 


0.012743(91) 


0.029966(79) 0.000674(77) 


6 


0.001423(76) 


0.005091(71) 


0.014570(57) 0.000136(77) 


7 


0.000469(53) 


0.002000(91) 


0.007103(46) 


8 


0.000155(43) 


0.000813(94) 


0.003505(69) 



Table 2: Average values of / x J Wilson loops of the actions A, M and B. 



action 


m 


ox 


c 


X 2 /d.o.f 


a(fm) 


A 


0.11(29) 


1.29(12) 


0.27(16) 


1.37 


0.52(5) 


M 


0.35(22) 


0.885(85) 


0.33(14) 


0.42 


0.43(4) 


B 


0.39(7) 


0.587(27) 


0.251(40) 


1.17 


0.35(2) 



Table 3: Results of the fits. The lattice spacing a is obtained by using ol = era? and y/cr = 427MeV. 
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Figure 1: Location of the actions in the (3u — (3\2 plane. Points A and B were obtained in [9]. 
Point M is interpolated at the mid-point between A and B. Symanzik and Iwasaki actions are also 
indicated by solid lines in the figure. 
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Figure 2: Comparison of 1 x T Wilson loop values between action A and the original blocked 
lattice. 
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Figure 3: Same as in Fig.2 but for action B. 
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Figure 4: Heavy quark potential for action A and the original lattice blocked at (3 = 6.00. The 
circles (diamonds) indicate on(off)-axis potentials. The curve in the figure is obtained by a fit to 
eq.(7). 
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Figure 6: Heavy quark potentials for actions A, B and M. The 3 lines almost on top of each other 
are separate fits for each action. 
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